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Abstract. The recent 23-30 January and 7-11 March 2012 
solar proton event (SPE) periods were substantial and caused 
significant impacts on the middle atmosphere. These were 
the two largest SPE periods of solar cycle 24 so far. The 
highly energetic solar protons produced considerable ion- 
ization of the neutral atmosphere as well as HO x (H, OH, 
HO 2 ) and NO x (N, NO, NO 2 ). We compute a NO x produc- 
tion of 1.9 and 2.1 Gigamoles due to these SPE periods in 
January and March 2012, respectively, which places these 
SPE periods among the 12 largest in the past 50 yr. Aura Mi- 
crowave Limb Sounder (MLS) observations of the peroxy 
radical, HO 2 , show significant enhancements of >0.9ppbv 
in the northern polar mesosphere as a result of these SPE pe- 
riods. Both MLS measurements and Goddard Space Flight 
Center (GSFC) two-dimensional (2-D) model predictions in- 
dicated middle mesospheric ozone decreases of > 20 % for 
several days in the northern polar region with maximum de- 
pletions >60% over 1-2 days as a result of the HO x pro- 
duced in both the January and March 2012 SPE periods. 
The SCISAT-1 Atmospheric Chemistry Experiment Fourier 
Transform Spectrometer (ACE) and the Envisat Michelson 
Interferometer for Passive Atmospheric Sounding (MIPAS) 
instruments measured NO and NO 2 (~NO x ), which indi- 
cated enhancements of over 20 ppbv in most of the northern 
polar mesosphere for several days as a result of these SPE pe- 


riods. The GSFC 2-D model and the Global Modeling Initia- 
tive three-dimensional chemistry and transport model were 
used to predict the medium-term (~ months) influence and 
showed that the polar middle atmospheric ozone was most af- 
fected by these solar events in the Southern Hemisphere due 
to the increased downward motion in the fall and early win- 
ter. The downward transport moved the SPE-produced NO y 
to lower altitudes and led to predicted modest destruction of 
ozone (5-13%) in the upper stratosphere days to weeks af- 
ter the March 2012 event. Polar total ozone reductions were 
predicted to be a maximum of 1.5% in 2012 due to these 
SPEs. 


1 Introduction 

Solar flares erupted in January and March 2012 and sent 

large fluxes of charged particles towards the Earth. Much 

of the energy was carried by solar protons, which impacted 
the middle atmosphere (stratosphere and mesosphere) lead- 

ing to ionizations, dissociations, dissociative ionizations, and 
excitations. Such periods are generally known as solar proton 
events (SPEs) and cause proton-induced atmospheric interac- 
tions resulting in the production of HO x (H, OH, HO 2 ), and 
NO x (N, NO, NO 2 ) constituents either directly or through 
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a photochemical sequence in the polar middle atmosphere. 
Such HO x and NO x enhancements can then lead to a pertur- 
bation of ozone. 

The major influence of these eruptions caused enhanced 
solar proton flux levels above background at the Earth for 
several days in January and March 2012. Although SPEs and 
their middle atmospheric constituent influences have been 
studied before (e.g., Swider and Keneshea, 1973; Heath et 
al., 1977; Solomon et al., 1981; McPeters and Jackman, 
1985; Randall et al., 2001; Lopez-Puertas et al., 2005a; von 
Clarmann et al., 2005; Verronen et al., 2006, 2008, 2011a, 
b; Krivolutsky et al., 2006; Funke et al., 2011; Jackman et 
al., 1990, 2001, 2005, 2007, 2008, 2009, 2011; Sinnhuber 
et al., 2012), each new event offers a chance to investigate 
the SPE-caused atmospheric perturbation under somewhat 
different conditions. Also, such events have been shown to 
cause some longer lasting stratospheric changes (e.g., Ran- 
dall et al., 2001; Jackman et al., 2008, 2009) and need to be 
considered within the context of other ongoing stratospheric 
changes (WMO, 2011). 

The recent study of von Clarmann et al. (2013) showed 
several atmospheric constituent changes due to the January 
and March 2012 SPEs using Envisat Michelson Interfer- 
ometer for Passive Atmospheric Sounding (MIPAS) instru- 
ment measurements. These latest large events, which pro- 
vided very significant middle atmospheric impulses over 
short (~ days) periods of time, also offer new opportunities 
to test the general understanding of an atmospheric change 
caused by very substantial perturbations. 

For example, the January 2012 SPE period occurred in the 
Northern Hemisphere (NH) winter and Southern Hemisphere 
(SH) summer. The solar zenith angle (SZA) was very differ- 
ent in the NH and SH polar regions for this month of the 
year. Since the SZA dramatically impacts the background 
atmosphere upon which the SPE-caused perturbation inter- 
acts, the ultimate SPE-caused atmospheric variation is highly 
dependent on the SZA. The influences of the SPE-produced 
HO x and NO x species impacted ozone somewhat differently 
in the two hemispheres in January. Also, the NH was subject 
to a sudden stratospheric warming (SSW) event along with 
a mesospheric cooling in January (e.g., von Clarmann et al., 
2013 and Paivarinta et al., 2013). 

The March 2012 SPE period occurred closer to equinox 
for the two hemispheres, however, the NH was transitioning 
to spring and the SH was transitioning to fall. This meant that 
the SPE-caused NO x was transported somewhat differently 
in the two hemispheres and had a larger impact on ozone in 
the SH than in the NH. 

We will show the measured influence of these SPEs on 
HO 2 and ozone from the Aura Microwave Limb Sounder 
(MLS) instrument. The influence of these SPEs on the 
sum of NO and NO 2 (~NO x ) from the SC1SAT-1 Atmo- 
spheric Chemistry Experiment Fourier Transform Spectrom- 
eter (ACE-FTS) and the Envisat MIPAS instruments will also 
be shown. These observations of the atmospheric response to 


the January and March 2012 SPE periods will be compared 
with global model predictions. 

This paper is divided into seven primary sections, includ- 
ing the Introduction. The solar proton flux and ionization rate 
computation are discussed in Sect. 2 and SPE-induced pro- 
duction of HO x and NO x are discussed in Sect. 3. A descrip- 
tion of the Goddard Space Flight Center two-dimensional 
model and the Global Modeling Initiative three-dimensional 
chemistry transport model is given in Sect. 4. Results from 
the two models for short-term (days) constituent changes, 
with comparisons to measurements for these SPEs, are 
shown in Sect. 5 while medium-term (months) constituent 
changes caused by these SPEs are discussed in Sect. 6. The 
conclusions are presented in Sect. 7. 


2 Proton flux and ionization rate 

We use the solar proton flux (energies 1 to 300 Me V) pro- 
vided by the National Oceanic and Atmospheric Administra- 
tion (NOAA) Space Weather Prediction Center (SWPC) for 
the NOAA Geostationary Operational Environmental Satel- 
lites (GOES) (see http://www.swpc.noaa.gov/ftpmenu/lists/ 
particle.html). The GOES 13 data are considered to be the 
most reliable of the current GOES datasets for the proton 
fluxes depositing energy into polar latitudes and were used as 
the source of protons in several energy intervals for the very 
active time periods 23-30 January 2012 and 7-11 March 
2012. For example, the GOES proton flux during the 22-30 
January 2012 time period is given in Fig. 1 (Top) for sev- 
eral different energies (> 1, > 10, > 30, > 50, and >100 MeV). 
Note that the proton flux levels increased by a factor of a 
thousand or more for energies > 10, >30, and >50 MeV on 
23 January 2012. Four days later (27 January) another event 
led to proton flux levels increasing by a factor of more than a 
hundred for > 30, > 50, and > 100 MeV. 

The proton flux data were used to compute the ion pair 
production profiles employing the energy deposition method- 
ology discussed in Jackman et al. (1980), where the creation 
of one ion pair was assumed to require 35 eV (Porter et al., 
1976). The SPE-produced daily average ionization rates for 
the polar cap regions (> 60° geomagnetic latitude) are given 
in Fig. 1 (Bottom) for a nine day period in January 2012 from 
lOhPa (~30km) to 0.001 hPa (~95km). Peak ionization 
rates above 2000 cm -3 s -1 on 24 January were computed for 
the 0.01 to 0.2 hPa region. The ionization was greatly dimin- 
ished by 26 January, but then the second SPE started late 
on 27 January and reached a peak on 28 January (also, see 
Fig. 1, Top). 

The GOES proton flux during the 6-11 March 2012 time 
period is given in Fig. 2 (Top) for the same proton ener- 
gies as in Fig. 1 (Top). The proton flux levels increased by 
a factor of a thousand or more for several energies (from 
>10 through >100 MeV) on 7-8 March 2012. The SPE- 
produced daily average ionization rates during this event for 
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GOES Proton Meas.: Jan. 2012 SPEs 



Ion rates (#cm 3 s '): Jan. 2012 SPEs 



January 201 2 


Fig. 1. (Top) GOES 13 proton flux measurements in January 2012 
for energies > 1 MeV (black), >10MeV (red), >30MeV (green), 
>50MeV (gold), and >100 MeV (blue). These data are provided 
by the NOAA SWPC at their website (see Sect. 2). ( Bottom ) Daily 
averaged ionization rates over the 10 to 0.001 hPa pressure range 
for the 22-30 January 2012 time period. Contour intervals are 100, 
200, 500, 1000, and 2000 (#cm -3 s -1 ). 

the polar cap regions are given in Fig. 2 (Bottom). Peak ion- 
ization rates above 2000 cm -3 s -1 on 8 March were com- 
puted for the 0.01 to 0.3 hPa region. Large ionization rates 
>100 cm -3 s -1 are computed all the way down to lOhPa 
on 7-8 March. As illustrated in Fig. 2 (Top) compared with 
Fig. 1 (Top), this SPE period consisted of a larger flux of pro- 
tons at higher energies (> 100 MeV) than were present in the 
January 2012 SPE period. The ionization is greatly dimin- 
ished by 1 1 March. 

3 HO x (H, OH, H0 2 ) and NO x (N, NO, NO z ) 

production 

Charged particle precipitation results in the production of 
HO x through complex positive ion chemistry (Solomon et al., 
1981). The charged particle-produced HO x is a function of 
ion pair production and altitude and is included in model sim- 
ulations using a lookup table from Jackman et al. (2005, Ta- 
ble 1), which is based on the work of Solomon et al. (1981). 
Each ion pair results in the production of about two HO x 
constituents for the upper stratosphere and lower mesosphere 
and less than two HO x constituents for the middle and upper 
mesosphere. 


GOES Proion Meas.: Mar. 2012 SPEs 




Fig. 2. (Top) GOES 13 proton flux measurements in March 2012 
for energies > 1 MeV (black), >10 MeV (red), >30 MeV (green), 
>50 MeV (gold), and >100 MeV (blue). These data are provided 
by the NOAA SWPC at their website (see Sect. 2). (Bottom) Daily 
averaged ionization rates over the 10 to 0.001 hPa pressure range for 
the 6-1 1 March 2012 time period. Contour intervals are 100, 200, 
500, 1000, and 2000 (#cm -3 s -1 ). 

NO x is produced when the energetic charged particles 
(protons and associated secondary electrons) dissociate N 2 
as they precipitate into the atmosphere. Here it is assumed 
that ~ 1 .25 N atoms are produced per ion pair and the proton 
impact of N atom production is divided between the ground 
state N( 4 S) (~ 45 % or ~ 0.55 per ion pair) and excited state 
N( 2 D) (~ 55 % or ~ 0.7 per ion pair) nitrogen atoms (Porter 
et al., 1976). We calculate that the 23-30 January 2012 SPE 
period produced about 1.9 Gigamoles of NO x and that the 
7-1 1 March 2012 SPE period produced about 2.1 Gigamoles 
of NO x . Thus, the January 2012 and the March 2012 SPE 
periods were computed to be the twelfth and tenth largest in 
the past 50 yr, respectively (see Jackman et al., 2008). 

4 Model predictions 

4.1 Description of the GSFC 2-D model 

The latest version of the Goddard Space Flight Center 
(GSFC) two-dimensional (2-D) atmospheric model was used 
to predict the impact of solar protons on the atmosphere. This 
model was first discussed over 20 yr ago (Douglass et al., 
1989; Jackman et al., 1990) and has undergone extensive 
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improvements over the years (e.g., Considine et al., 1994; 
Jackman et al., 1996; Fleming et al., 1999, 2007, 2011). The 
vertical range of the model, equally spaced in log pressure, is 
from the ground to approximately 92 km (0.0024 hPa) with a 
1 km grid spacing. The model has a 4° latitude grid spacing. 

For this study, the transport is computed off-line and is de- 
rived using the daily average global winds and temperatures 
from the NASA Modern Era Retrospective-analysis for Re- 
search and Applications (MERRA) meteorological analysis 
(see the website: http://gmao.gsfc.nasa.gov/researcli/merra/) 
for 1979-2012. Thirty-day running averages of the residual 
circulation, eddy diffusion, zonal mean wind, and zonal mean 
temperature are computed using the methodology detailed in 
Fleming et al. (2007), and are used as input into the GSFC 
2-D model. For this paper we used the transport fields for the 
year 2012, which included the January SSW. 

The ground boundary conditions in the GSFC 2-D model 
for the source gases are taken from WMO (2011) for year 
2012. The model uses a chemical solver described in Jack- 
man et al. (2005) and Fleming et al. (2007, 2011). The photo- 
chemical gas and heterogeneous reaction rates and photolysis 
cross sections have been updated to the latest Jet Propulsion 
Laboratory recommendations (Sander et al., 2010) for these 
computations. 

4.2 Simulations of the GSFC 2-D model 

The starting conditions for the two GSFC 2-D model sim- 
ulations used in this study were provided by a time de- 
pendent “spin-up” simulation from January 2000-December 
2011, which included the appropriate source gas boundary 
conditions from WMO (2011) and the changing transport 
fields derived from MERRA for this time period. Two time- 
dependent simulations, “A” and “B”, were continued over the 
period 1 January through 31 December 2012. Simulation “A” 
included no SPEs and simulation “B” included all SPEs in 
this period. 

4.3 Description of the GMI 3-D CTM 

The Global Modeling Initiative (GMI) three-dimensional (3- 
D) chemistry and transport model (CTM) was also used in 
this study. The GMI 3-D CTM is well suited to simulate the 
SPE-produced NO x in the lower mesosphere and through- 
out the stratosphere during 2012. The 1° x 1.25° horizontal 
resolution of the latest version of the GMI 3-D CTM (Stra- 
han et al., 2013) uses MERRA meteorological fields for the 
time periods under investigation. More information about the 
GMI 3-D CTM and its chemical mechanism can be found in 
Strahan et al. (2007) and Duncan et al. (2007). The strato- 
spheric chemical mechanism is taken from Douglass and 
Kawa (1999) and the polar stratospheric cloud parameteri- 
zation is described in Considine et al. (2000). 

The vertical range of the GMI 3-D CTM is from the 
ground up to approximately 80 km (0.015 hPa) and the model 


lacks some mesospheric processes, such as Lyman alpha 
photochemistry. Thus, the GMI 3-D CTM cannot be used to 
address the short-term influence of SPE-produced HO x and 
its impact on ozone in the middle and upper mesosphere. The 
GMI 3-D CTM will, therefore, be used in conjunction with 
the GSFC 2-D model to study the medium-term influence of 
the 2012 SPEs. 

4.4 Simulations of the GMI 3-D CTM 

The starting conditions for the two GMI 3-D CTM sim- 
ulations used in this study were provided by a time de- 
pendent “spin-up” simulation from January 2004-December 
2011, which included the appropriate source gas boundary 
conditions from WMO (2011) and the changing transport 
fields derived from MERRA for this time period. Two time- 
dependent simulations, “C” and “D”, were continued over 
the period 1 January through 31 December 2012. Simulation 
“C” included no SPEs and simulation “D” included all SPEs 
in this period. 

5 Short-term Influences of the SPEs 

A number of constituents have been observed and modeled 
to be modified by SPEs (e.g., Lopez-Puertas et al., 2005a, 
b; von Clarmann et al., 2005; Verronen et al., 2006; Jack- 
man et al., 2011; Damiani et al., 2012). We will focus on the 
SPE-caused impact on HO x , NO x , and ozone in this study. 
The SPEs can cause substantial increases in HO x and NO x 
constituents (see Sect. 3), which can directly impact ozone 
through catalytic destruction cycles (e.g., see Johnston and 
Podolske, 1978). 

5.1 Peroxy radical (HO 2 ) 

The Aura Microwave Limb Sounder (MLS) measurement of 
the peroxy radical (HO 2 ) provides a direct observation and 
quantification ofHO x changes due to the January 2012 SPE 
in the middle atmosphere. Figure 3 (Top left) shows the daily 
average MLS HO 2 changes from a five-day (18-22 January) 
average of HO 2 measurements before the SPE period for the 
60-82.5° N band. Enhancements of >0.2ppbv are observed 
in the middle mesosphere over periods 23-26 January and 
28-29 January. A peak of > 1 ppbv is observed from ~ 0.05- 
0.1 hPa on 24 January. These HO 2 enhancements are slightly 
larger than those reported in Jackman et al. (2011) for the 
January 2005 SPEs. 

The measurement of HO 2 in the atmosphere by MLS is 
impressive given its relatively small abundance throughout 
the middle atmosphere, however, some care must be exer- 
cised in using these data. For example, MLS HO 2 has a bias 
that is usually removed by taking day-night differences over 
the entire recommended pressure range. It was advantageous 
in this work to include all the MLS HO 2 measurements to 
determine the SPE-caused impact. Removal of the five-day 
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average before the SPE period, as done in this study, will 
also lead to a removal of the MLS HCb bias assuming that 
the bias varies in only a minor way over a couple of weeks. 
Although the MLS HO 2 measurements are not recommended 
for scientific use above 0.046 hPa, those MLS observations 
are shown here to qualitatively indicate the morphology of 
the HO 2 change during the disturbed periods. 

We also show the daily average MLS HO 2 changes from 
a five-day (18-22 January) average of HO 2 measurements 
before the SPE period for the 60-82.5° S band in Fig. 3 (Top 
right). Note that the observed SH HO 2 enhancement is at a 
peak value of >0.5 ppbv for the ~ 0.015-0.07 hPa region on 
24 January, which is less than the peak levels shown in the 
NH (see Fig. 3, Top left). The HO x produced by the SPE is 
approximately the same in both hemispheres, however, the 
lifetime of HO 2 is diminished in the SH (summer) due to 
the higher background HO x level and a faster loss of HO 2 
at this time of year. The increased solar radiation in the SH 
(summer) leads to an increase in the amount of background 
mesospheric HO x by a factor of about 2-6 relative to the NH, 
thus the loss of HO 2 is larger in the SH through the increase 
in the rate of reactions 

0H + H0 2 ->H 2 0 + 0 2 (Rl) 

and 

H + HO 2 — »■ products (R2) 

The GSFC 2-D model predicted HO 2 change from the five- 
day (18-22 January) average using simulation “B” is shown 
in Fig. 3 (Bottom plots) for the same period. Also, the MLS 
averaging kernel (AK) is applied to the model results. The 
variation with time and pressure as well as the magnitude of 
the computed change is fairly similar to the MLS HO 2 obser- 
vations, especially in the NH. The predicted HO 2 increase is 
slightly smaller than the observed increase due to the SPEs in 
the SH, however, the pressure levels impacted are fairly sim- 
ilar between the model and measured values. The predicted 
HO 2 enhancements due to the SPEs are within the computed 
precision error (~ 0.4 ppbv near 0.1 hPa) for daily average 
MLS measurements of HO 2 in these latitude bands. 

We also completed a similar analysis of the HO 2 enhance- 
ments caused by the March 2012 SPEs. The HO 2 changes 
from a five-day (2-6 March) average of HO 2 measurements 
before the SPE period for the 60-82.5° N and 60-82.5° S 
bands were examined. We found MLS observed enhance- 
ments of >0.2 ppbv in the mesosphere over five days (7-11 
March) with a peak of 0.8 ppbv occurring near 0.05 hPa on 8 
March in the NH and observed enhancements >0.2 ppbv on 
three days (7, 8, and 10 March), see Fig. 4 (Top left), and a 
peak of 0.4 ppbv occurring near 0.05 hPa on 8 March in the 
SH (Fig. 4, Top right). 

The general location of the peaks and variations with pres- 
sure are fairly similar between the measurements and model 
predictions for HO 2 in early March 2012 during the SPE, 


H0 2 ppbv change in Jan. 2012 



2D-AK (60-82. 5°N) 2D-AK (60-82. 5°S) 



January 2012 January 2012 


Fig. 3. Daily averaged HO 2 changes from Aura MLS measurements 
{Top) and GSFC 2-D model predictions {Bottom) for the 60-82.5° N 
band {left plots) and 60-82.5° S band {right plots). An average ob- 
served (predicted) HO 2 profile for the period 18-22 January 2012 
was subtracted from the observed (predicted) HO 2 values for the 
plotted days of 22-30 January 2012. The contour intervals for the 
HO 2 differences are —0.2, —0.1, 0.0, 0.1, 0.2, 0.5, and 1 ppbv. The 
MLS averaging kernel (AK) was used to sample the model results. 

however, there are some differences in the peak values (see 
Fig. 4, Bottom plots). The predicted HO 2 change as a result 
of the March 2012 SPE shows a peak on 8 March near 0.02 to 
0.05 hPa of about 1 ppbv for latitude band 60-82.5° N and a 
peak of about 0.7 ppbv for latitude band 60-82.5° S. The so- 
lar radiation is slightly larger in the SH compared to the NH 
polar region in early March. Also, more mesospheric H 2 O is 
present in the SH compared to the NH due to the upward mo- 
tion during the southern polar summer transporting up H 2 O 
from the upper stratospheric water source. These two fac- 
tors lead to a larger background concentration of HO x con- 
stituents in the SH and thus a larger loss of the SPE-produced 
HO 2 through Rl and R2 and a slightly smaller peak of HO 2 
in the SH compared to the NH. 

5.2 Ozone 

The HO x constituents are the primary cause of ozone deple- 
tion in the mesosphere during SPEs and the impact during 
the January and March 2012 SPE periods is seen clearly in 
the MLS ozone measurements. Several catalytic HO x cycles 
are important in the middle atmosphere, however, the process 

H + 0 3 -*0H + 0 2 (R3) 
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H0 2 ppbv change in Mar. 2012 
MLS (60-82. 5°N) MLS (60-82. 5°S) 



80 
70, 
60 ; 
50 


30 


2D-AK (60-82. 5°N) 


2D-AK (60-82. 5°S) 



March 201 2 


March 2012 


Fig. 4. Daily averaged H0 2 changes from Aura MLS measurements 
{Top) and GSFC 2-D model predictions {Bottom) for the 60-82.5° N 
band {left plots) and 60-82.5° S band {right plots). An average ob- 
served (predicted) HO 2 profile for the period 2-6 March 2012 was 
subtracted from the observed (predicted) HO 2 values for the plotted 
days of 6-1 1 March 2012. The contour intervals for the HO 2 differ- 
ences are —0.5, —0.2, —0.1, 0.0, 0.1, 0.2, 0.5, and 1 ppbv. The MLS 
averaging kernel (AK) was used to sample the model results. 


Ozone % change 
MLS ( 60-82. 5°N) 


n Jan. 201 2 

MLS (60-82. 5°S) 



22 23 24 25 26 27 28 29 30 31 
January 2012 


22 23 24 25 26 27 28 29 30 31 
January 2012 


Fig. 5. Daily averaged ozone changes from Aura MLS measure- 
ments {Top) and GSFC 2-D model predictions {Middle, Bottom) 
for the 60-82.5° N band {left plots) and the 60-82.5° S band {right 
plots). An average observed (predicted) ozone profile for the period 
18-22 January 2012 was subtracted from the observed (predicted) 
ozone values for the plotted days of 22-30 January 2012 for Aura 
MLS measurements {Top) and for the GSFC 2-D model simulation 
"B (with SPEs)” {Middle). The Bottom plot shows the difference be- 
tween GSFC 2-D model simulations “B (with SPEs)” and “A (with- 
out SPEs).” The contour intervals for the ozone differences are —80, 
-60, -40, -20, -10, -5, -2,-1, 0, 1, 2, 5, 10, 20, and 40%. The 
MLS averaging kernel (AK) was used to sample the model results. 


followed by 

0H + 0-*H + 0 2 (R4) 

Net: 0 + 0 3 -* 20 2 (R5) 

is dominant throughout much of the mesosphere. During 
strong SPEs, such as these in 2012, the HO x catalytic cycle 

OH + 0 3 -* H0 2 + 0 2 (R6) 

followed by 

H0 2 + 0-> OH + 0 2 (R7) 

Net: O + O3 -* 20 2 (R8) 

also contributes significantly to the ozone decrease. 

We show the MLS daily average ozone change from a 


five-day (1 8-22 January) average of ozone measurements be- 
fore the SPE period for the 60-82.5° N band in Fig. 5 (Top 
left). Ozone depletions of > 20 % are observed in the middle 
mesosphere for the 60-82.5° N band over the 24-27 January 
period with a maximum of 60% depletion occurring near 
0.02-0.05 hPa on 25 January (Fig. 5, Top left). The HO x con- 
stituents have a relatively short lifetime (~ hours) throughout 
most of the mesosphere (below ~ 80 km, see Pickett et al., 
2006), however, they can cause very large ozone depletion. 


which has also been observed in other very large substantial 
SPEs (e.g., Solomon et al., 1983; Jackman et al., 2001, 2011; 
Verronen et al., 2006). 

Surprisingly, the MLS observations show an increase in 
ozone at their retrieval level of 0.046 hPa for 28-30 January, 
with ozone decreases above and below. It is unclear what 
causes this MLS-observed ozone increase at the 0.046 hPa 
level during this period. M1PAS ozone measurements at this 
same time period and location also show this behavior (von 
Clarmann et al., 2013). The observed ozone increases at pres- 
sures greater than 0.5 hPa are caused by seasonal changes on- 
going at this time of year, which mitigate and overwhelm the 
SPE-caused ozone decrease. 

We show two model computations of the SPE-caused pre- 
dicted ozone change in Fig. 5 (Middle and Bottom left). The 
model predicted changes in Fig. 5 (Middle left) are derived in 
the same manner as the MLS observed changes, wherein the 
modeled ozone in simulation “B” is subtracted from the five- 
day (18-22 January) average of modeled ozone before the 
SPE period. Thus, this model predicted change will also have 
seasonal variations included. The model predicted changes 
for Fig. 5 (Bottom left) are derived by subtracting simula- 
tion “A” from simulation ”B” to produce a SPE-only caused 
ozone change thereby removing the seasonal ozone increase 
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in the upper stratosphere. As in Figs. 3 and 4, the MLS aver- 
aging kernel (AK) is applied to the model results. Note that 
there are some differences between the measurements and 
model predictions (see Fig. 5, Top and Middle left). For ex- 
ample, although both the measurements and model predic- 
tions indicate a SPE-caused ozone depletion in the middle 
mesosphere over the 24-27 January period, the model pre- 
dicts a somewhat larger altitude range of predicted ozone de- 
crease from the SPEs as well as a smaller seasonal increase 
in the upper stratosphere and lower mesosphere. 

As in Fig. 5 (Top left), we show the MLS daily aver- 
age ozone change from a five-day (18-22 January) aver- 
age of ozone measurements before the SPE period for the 
60-82.5° S band in Fig. 5 (Top right). Ozone depletions of 
>10% are observed in the middle mesosphere for the 60- 
82.5° S band over the 23-24 January period with a maximum 
of 30% depletion occurring near 0.02 hPa on 23 January. A 
small depletion from about 2-0.04 hPa is even observed on 
28 January as a result of the second SPE in January. These 
SPE-caused ozone depletions are less in the polar SH (com- 
pared to the polar NH, left plots) due to the shorter lifetime 
of HO x (and ultimate HO x influence) in the summer relative 
to the winter (e.g., see Solomon et al., 1983; Jackman et al., 
2008). 

We show two computations of the SPE-caused predicted 
ozone change for the 60-82.5° S band in Fig. 5 (Middle and 
Bottom right), derived in the same way as the model predic- 
tions in Fig. 5 (Middle and Bottom left). The model predicted 
ozone change is in reasonable agreement with the observed 
MLS ozone changes on 23-24 January, both in the magnitude 
and timing. The model predicted changes shown in Fig. 5 
(Bottom right), derived by subtracting simulation “A” from 
“B”, indicate that seasonal increases in ozone are important 
in the upper stratosphere (~ 1-6 hPa) and above 0.2 hPa on 
29-31 January. 

Given the somewhat surprising MLS observations of an 
increase in ozone at the retrieval level of 0.046 hPa in the 60- 
82.5° N band for the SPE on 28-30 January (see Fig. 5, Top 
left), we show the MLS observations during the March 2012 
SPE period for this same region in Fig. 6 (Top left). Clearly, 
MLS measured mesospheric ozone is being depleted at all 
levels as a result of the SPE. The model results given in Fig. 6 
(Middle and Bottom left), computed in the same way as in 
Fig. 5, also show ozone being depleted throughout the meso- 
sphere. Both measurements and model results show ozone 
decreases > 60 % on 8-9 March at about 0.05-0.02 hPa. The 
measurements and model simulations also indicate a modest 
seasonally-driven ozone increase in the upper stratosphere. 
As an aside, MIPAS observations show similar ozone deple- 
tion due to the March 2012 SPE. 

We also examined the MLS ozone observations as well 
as our model predictions for the 60-82.5° S band during the 
March 2012 SPE period (see Fig. 6, Top right). The mea- 
surements and model results indicate significant ozone de- 
pletion as a result of the SPE with both showing ozone de- 
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Fig. 6. Daily averaged ozone changes from Aura MLS measure- 
ments (Top) and GSFC 2-D model predictions ( Middle , Bottom) for 
the 60-82.5° N band ( left plots) and for the 60-82.5° S band ( right 
plots). An average observed (predicted) ozone profile for the pe- 
riod 2-6 March 2012 was subtracted from the observed (predicted) 
ozone values for the plotted days of 6-11 March 2012 for Aura 
MLS measurements (Top) and for the GSFC 2-D model simulation 
"B (with SPEs)” (Middle). The Bottom plot shows the difference be- 
tween GSFC 2-D model simulations ”B (with SPEs)” and ‘A (with- 
out SPEs).” The contour intervals for the ozone differences are —80, 
-60, -40, -20, -10, -5, -2, -1, 0, 1, 2, 5, and 10%. The MLS 
averaging kernel (AK) was used to sample the model results. 

creases > 20 % on 8 March through the middle mesosphere. 
The model predictions given in Fig. 6 (Middle and Bottom 
right) do indicate a somewhat larger ozone depletion for the 
7-11 March period, when compared to the MLS measure- 
ments. 

5.3 Nitrogen oxides, NO x (NO+NO 2 ) 

ACE-FTS (hereinafter referred to as ACE) (Bemath et al., 
2005) and MIPAS (Fischer et al., 2008) measured both NO 
and NO 2 (e.g., Rinsland et al., 2005; Lopez-Puertas et al., 
2005a). We defined NO x in Sect. 1 as N + NO + NO 2 , but 
here we identify NO x as approximated by NO + NO 2 since 
N is minuscule in the middle atmosphere. 

ACE provided NO + NO 2 profile measurements at fairly 
high latitudes during the January through April 2012 period. 
For example, for the period Days of Year (DoY) 22-47 (22 
January to 16 February) ACE observations were made in the 
NH between about 65 and 68° N. The ACE measured NO x 
changes from DoY 22 (22 January) are shown in Fig. 7a. 
Large increases (>100ppbv) are observed in the middle to 
upper mesosphere (<0.05hPa) for DoY 24-36. Some of 
these increases are certainly due to the solar protons and 
associated secondary electrons, however, energetic electron 
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Fig. 7. Daily averaged NO x (NO + NO 2 ) changes from Day of 
Year (DoY) 22 (22 January) 2012 values as measured by ACE (a) 
and MIPAS (65-68° N) (b). The contour intervals for the NO x dif- 
ferences are 1, 2, 5, 10, 20, 50, 100, 200, and 500 ppbv. The northern 
polar latitudes of ACE measurements are given in the Top plot. 


precipitation (EEP) may also be contributing in the mid- 
dle and upper mesosphere (pressures <0.1 hPa). The solar 
protons penetrate deeper into the atmosphere than the EEP, 
and cause most of the measured enhanced NO x in the lower 
mesosphere. 

Substantial increases (> 20 ppbv) were observed through- 
out most of the mesosphere for these same days continuing 
through Day of Year 44 for the lower mesosphere (~ 0.04- 
0.7 hPa). MIPAS measurements are also available during this 
time period, but not for every day, and are given in Fig. 7b. 
It should be noted that the MIPAS nominal observations 
used here are limited to tangent heights of 70 km. Hence, 
the retrieved NO x profiles have little-to-no vertical resolu- 
tion above 70 km. In the SH, MIPAS NO x might be underes- 
timated (due to the effect of the a priori profile shape) even 
down to 65 km. 

In spite of these important caveats, the MIPAS measure- 
ments are generally in agreement with ACE when the slightly 
larger averaging kernels of MIPAS are taken into account. 
Such measured NO x enhancements were similar to those ob- 
served by ACE as a result of the January 2005 SPEs (Jack- 
man et al., 2011). 

ACE also takes measurements at high southern latitudes, 
however, the number of days of observations is relatively 
limited and will not be shown. A comparison of the SPE- 
caused NO x enhancement over a larger polar region (60-90°) 
is shown in Fig. 8 for the NH (Top left) and SH (Top right), 
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Fig. 8. Daily averaged NO x (NO + NO 2 ) changes from DoY 22 
(22 January) 2012 values for 60-90° N ( left plots) and 60-90° S 
(right plots) as measured by MIPAS (Top) and calculated by the 
GSFC 2-D model (Middle) and the GMI 3-D CTM (Bottom). The 
contour intervals for the NO x differences are 1, 2, 5, 10, 20, 50, and 
100 ppbv. 


respectively, for the MIPAS instrument and the GSFC 2-D 
model (Middle plots) and the GMI 3-D CTM (Bottom plots). 
The GMI 3-D CTM results are most reliable at pressures 
greater than 0.04 hPa, thus all plots in Fig. 8 are only shown 
up to this level. The models indicate a somewhat smaller NO x 
change in the NH and a somewhat larger NO x change in the 
SH, compared with MIPAS observations. The smaller pre- 
dicted NO x change in the NH may be related to the absence 
of concurrent energetic electron precipitation in our compu- 
tations. 

As an aside, it is difficult to include the EEP contribu- 
tion in our computations during the time period of the SPEs. 
The NOAA MEPED instrument’s electron energy channels 
are contaminated by protons, thus any measurements of EEP 
during SPEs is questionable (e.g., p. 2 of Verronen et al., 
2011a; Table 3 of Yando et al., 2011). 

Although similar, there are some differences between the 
two model predictions of the enhanced NO x evolution. The 
GMI 3-D CTM holds the SPE-caused NO x increases closer 
to the NH pole (Fig. 8, Bottom left) and transports a larger 
amount of NO x to the stratosphere than indicated in the 
GSFC 2-D model (Fig. 8, Middle left) or observed in MI- 
PAS data (Fig. 8, Top left). In contrast, the GSFC 2-D model 
appears to conserve NO x over a longer period of time in the 
SH polar region than indicated in the GMI 3-D CTM (Fig. 8, 
Middle and Bottom right). 

We turn our attention to the March 2012 SPE period in 
Fig. 9. The ACE and MIPAS measured NO x ppbv change 
from 6 March is shown in Fig. 9a, b and indicates enhance- 
ments of over 20 ppbv throughout most of the northern polar 
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(a) NOx ppbv change (NH) — ACE 




Fig. 9. Daily averaged NO x (NO + NO 2 ) changes from DoY 66 
(6 March) 2012 values as measured by ACE (a) and MIPAS (78- 
83° N) (b). The contour intervals for the NO x differences are 1, 2, 
5, 10, 20, 50, and 100 ppbv. The northern polar latitudes of ACE 
measurements are given in the Top plot. 


mesosphere for several days (DoY 68-74 for ACE and DoY 
68-72 for MIPAS) as a result of this SPE period. Peak NO x 
increases over 100 ppbv were observed by ACE on DoY 68 
(8 March) between 0.01 and 0.02 hPa as a result of the March 
2012 SPE. 

As during the January 2012 SPE period, ACE also takes 
measurements at high southern latitudes. Again, however, 
the number of days of observations is relatively limited and 
will not be shown. A comparison of the SPE-caused NO x 
enhancement over a larger polar region (60-90°) is shown in 
Fig. 10 for the NH (Top left) and SH (Top right), respectively, 
for the MIPAS instrument and the GSFC 2-D model (Middle 
plots) and the GMI 3-D CTM (Bottom plots). Although the 
impacted altitudes are similar between the measurements and 
the model predictions, the models indicate a somewhat larger 
NO x change in both hemispheres, when compared with MI- 
PAS observations. Some of this larger model predicted SPE- 
caused NO y enhancement can be explained by the finer ver- 
tical grid of the model ( 1-2 km in the mesosphere) compared 
to the MIPAS vertical field-of-view of about 3 km (Lopez- 
Puertas et al., 2005a). As for January (see Fig. 8 ), the GMI 
3-D CTM transports a larger amount of NO x to the polar 
NH stratosphere (Fig. 10, Bottom left) than indicated in the 
GSFC 2-D model (Fig. 10, Middle left). The models show 
generally similar NO x transport to the stratosphere in the SH 
(Middle and Bottom right). 

Both models predict that the initial NO x enhancement is 
substantially larger than the measured enhancements, with 
peak increases greater than 50 ppbv on DoY 68-69 in the up- 
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Fig. 10. Daily averaged NO x (NO + NO 2 ) changes from DoY 66 
(6 March) 2012 values for 60-90° N (left plots) and 60-90° S (right 
plots) as measured by MIPAS (Top) and calculated by the GSFC 
2-D model (Middle) and the GMI 3-D CTM (Bottom). The contour 
intervals for the NO x differences are 1, 2, 5, 10, 20, and 50 ppbv. 


per mesosphere for both hemispheres, compared with peak 
MIPAS measured increases of only about 25 ppbv (Fig. 10, 
left). Note also that both models predicted NO x changes in 
the NH on DoY 76-80 showing >5 ppbv throughout most 
of the mesosphere (pressures <0.5hPa) whereas the MI- 
PAS measured NO x changes are < 5 ppbv in the same region 
(Fig. 10, left). MIPAS measured mesospheric NO x enhance- 
ments are gone by DoY 81, however, model predicted NO x 
increases continue and gradually diminish in the lower meso- 
sphere through DoY 97. The disappearance of NH MIPAS 
NO x enhancements is possibly related to a polar vortex rup- 
ture and release of mesospheric NO x towards mid latitudes, 
which appears in MIPAS CH 4 data (not shown). Both mod- 
els used transport from MERRA for the year 2012. However, 
since observational information in the mesosphere is limited 
there may be dynamical differences between the atmosphere 
and the models for the mesosphere. The models likely under- 
estimate the amount of high latitude mesospheric NO x trans- 
ported to lower latitudes. 

The model predicted NO x changes in the SH on DoY 76- 
80 are ~ 5-20 ppbv throughout the middle to lower meso- 
sphere (~ 0.03-0.7 hPa), which are in reasonable agree- 
ment with the MIPAS measurements during this time period 
(Fig. 10, right). In the upper stratosphere (~ 1-3 hPa) some- 
what more SPE-produced NO x is simulated in the models 
than observed. Thus, the models predict a more substantial 
SPE-caused NO x perturbation in the SH, than indicated in 
the measurements, which will then cause a larger impact on 
stratospheric ozone. 

Overall both models predict NO x enhancements which 
are larger than the MIPAS observations in both hemispheres 
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for the March 2012 SPE. Neither model includes ion chem- 
istry, thus both models will predict a smaller enhancement 
of HNO 3 (and subsequent reduction in NO x ) due to the 
SPEs. This repartitioning through ion chemistry of NO x to 
HNO 3 has been found to be an important process (see Lopez- 
Puertas et al., 2005b; Verronen et al., 2008, 2011b) and may 
be responsible for the loss of up to a few ppbv of NO x to 
HNO 3 in the lower mesosphere. This process is unlikely to 
explain the large differences between MIPAS and the models 
for NO x above ~ 0.1 hPa, but may reduce these discrepancies 
in the lower mesosphere. Also, the proton flux source for the 
March 2012 SPE is the same as that for the J anuary 20 1 2 SPE 
(GOES 13), thus this NO x over-prediction above ~0.1 hPa 
is not easily explained. Further model-measurement compar- 
isons of these SPE-caused impacts will be needed to help 
understand this disagreement. 


6 Medium-term Influences of the SPEs: odd nitrogen 
(NO y ) and ozone changes 

The last two figures (Figs. 11-12) show only GSFC 2- 
D model and GMI 3-D CTM predicted medium-term 
(~ months) influences of the SPE for all of 20 1 2. These plots 
were derived by subtracting the simulations without the SPEs 
(“A” for the GSFC 2-D model; “C” for the GMI 3-D CTM) 
from the simulations with the SPEs (“B” for the GSFC 2-D 
model and “D” for the GMI 3-D CTM) to produce a SPE- 
only caused impact. 

The NO x family rapidly converts in the stratosphere to 
other constituents in the odd nitrogen group (NO y = N( 4 S) 
+ N( 2 D) + NO + N0 2 + NO 3 + 2 N 2 O 5 + HNO 2 + HNO 3 + 
HO 2 NO 2 + CIONO 2 + BrON0 2 ). The SPE-caused NO y 
ppbv enhancement and ozone percentage changes are given 
in Figs. 1 1 and 12 for the polar (60-90°) NH and SH, respec- 
tively. NO y enhancements exceeding 2 ppbv are computed in 
the polar NH for much of the first five months of the year 
(J-F-M-A-M) throughout most of the mesosphere and part 
of the upper stratosphere in the GSFC 2-D model (Fig. 11, 
Top left). The GMI 3-D CTM predictes a greater impact by 
the 2012 SPEs on NO y in the polar NH (Fig. 11, Top right). 
Much more NO y is transported to the stratosphere from the 
mesosphere by the GMI 3-D CTM compared with the GSFC 
2-D model in the J-F-M time period. Some of that NO y en- 
hancement computed by the GMI 3-D CTM is diminished 
over the spring and summer, however, levels of NO y > 2 ppbv 
are observed for much of the year in the upper stratosphere 
(Fig. 11, Top right). 

Large computed NH ozone decreases (>20%, dark blue 
colors) during and shortly after the two SPE periods (Jan- 
uary and March) in the mesosphere were primarily caused 
by the SPE-produced HO x increases (Fig. 1 1, Bottom plots). 
The GSFC 2-D model predicted upper stratospheric ozone 
decreases of 2-4% (Fig. 11, Bottom left) and the GMI 3-D 
CTM predicted upper ozone decreases of 5-13% (Fig. 11, 
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Fig. 11. (Left) Difference between GSFC 2-D model simulation 
"B (with SPEs)” and “A (without SPEs)”. (Right) Difference be- 
tween GMI 3-D CTM simulation "D (with SPEs)” and “C (without 
SPEs)”. NO y (Top) and Ozone (Bottom) for 60-90° N are given for 
all months of 2012. The contour intervals for the NO y differences 
are 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, and 50 ppbv. The contour intervals 
for the Ozone differences are —20, —10, —5, —2, —1,0, 1, and 2 %. 


Bottom right). These differences in model predicted ozone 
depletions were caused by the somewhat larger transported 
NO y to the stratosphere given by the GMI 3-D CTM com- 
pared with the GSFC 2-D model (Fig. 11, Top plots). 

The primary catalytic cycle for NO y destruction of ozone 
is: 


NO + O 3 — > NO 2 + O 2 

(R9) 

followed by 


NO 2 + O — > NO + O 2 

(R10) 

Net : O + O 3 — > 2 O 2 . 

(Rll) 


There is a region of predicted enhanced ozone due to the 
SPEs in the middle to lower stratosphere (especially, see 
Fig. 11, Bottom right). The GMI 3-D CTM shows increased 
ozone of 1-2 % caused by the enhanced SPE-produced NO y 
sequestering chlorine and bromine into reservoir species 
(CIONO 2 and BrON 02 ). This reduces ozone depletion due 
to halogens in this region of the stratosphere. This mecha- 
nism ultimately leads to increased ozone in a portion of the 
middle and lower stratosphere and has been discussed before 
in relation to SPE impacts (e.g., see Jackman et al., 2009). 

The computed January-February NO y enhancements are 
not as substantial in the polar SH (Fig. 12) as in the polar 
NH due to considerable loss of odd nitrogen in the intense 
sunlight of SH summer and early fall. The reactions most 
important for this loss are: 

NO + /fu(< 191mn)->N + 0 (R12) 
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Fig. 12. (Left) Difference between GSFC 2-D model simulation 
“B (with SPEs)” and “A (without SPEs)”. (Right) Difference be- 
tween GMI 3-D CTM simulation "D (with SPEs)” and “C (without 
SPEs)”. NOy (Top) and Ozone (Bottom) for 60-90° S are given for 
all months of 2012. The contour intervals for the NOy differences 
are 0.1, 0.2, 0.5, 1, 2, 5, 10, and 20 ppbv. The contour intervals for 
the Ozone differences are —20, — 10, —5, —2, —1, and 0 %. 


followed by 

N + NO-*N 2 +0. (R13) 

The winds are also generally upwards in the SH upper strato- 
sphere in January and February, thus NO y constituents are 
transported upwards to the higher altitudes, where the loss 
process is greater. The vertical wind direction starts to reverse 
direction in early March and by late March is basically down- 
wards (Fig. 12, Top plots). The wind reversal drags much 
more NO y from the mesosphere to the upper stratosphere by 
late May (Fig. 12, Top plots) leading to NO y enhancements 
>2 ppbv and associated ozone depletion of 5-8 % (GSFC 2- 
D model) and 5-1 1 % (GMI 3-D CTM). NO y enhancements 
of > 0.2 ppbv and the associated ozone decreases of > 1 % de- 
scend through the winter and spring, reaching about 50hPa 
by November in both models. The predicted ozone depletion 
in the SFl can be thought of as an “upper limit” given that the 
predicted SPE-caused NO x increases were somewhat higher 
than the measured SPE-caused NO x increases (see Fig. 10, 
Right). Note, also, that there were other SPE periods in May, 
June, July, September, October, November, and December, 
which are seen in the enhanced NO y in both hemispheres 
(Figs. 11 and 12, Top panels). These particular SPE peri- 
ods were not as large as those that occurred in January and 
March. 

The impact of the 2012 SPEs on total ozone was computed 
in both models to be largest at the highest SH latitudes, where 
reductions of a maximum of about 1 % (GSFC 2-D model) to 


about 1.5% (GMI 3-D CTM) are computed several months 
after the March 2012 SPE period. The maximum total ozone 
changes are not predicted to occur during the SPE periods 
when the HO x catalytic cycles dominate; rather, the trans- 
port of the enhanced NO y to lower altitudes (and high am- 
bient ozone amounts) causes more total ozone impact. Total 
ozone reductions in the NH are predicted to be >0.5 % less 
than those in the SH. Thus, the impact of these SPE periods 
on total ozone is fairly small and it is unlikely that the per- 
turbation could be detected in total ozone satellite or ground 
measurements. 

Measurements show a 1-2 % total ozone increase from so- 
lar min to solar max due to solar ultraviolet radiation changes 
(e.g., Chandra and McPeters, 1994; Souldiarev and Hood, 
2006; Swartz et al., 2012) with fairly similar changes at low 
and high latitudes. The model computed high Southern po- 
lar latitude (70-90° S) impact from the 2012 SPEs on total 
ozone was found to be a decrease of >0.5 % for about seven 
months. The SPE-caused ozone depletion could thus be com- 
petitive with solar radiation-induced variations over periods 
of several months at polar latitudes. Also, the SPE impact is 
of opposite sign to the solar radiation effect: more SPEs oc- 
cur near solar max and lead to ozone decreases whereas so- 
lar radiation variation near solar max lead to ozone increases. 
SPE activity during a solar cycle therefore should be consid- 
ered carefully when examining the solar cycle induced vari- 
ations on polar ozone. 


7 Conclusions 

The January and March 2012 SPE periods caused substantial 
perturbations to the Earth’s middle atmosphere. In particular, 
measurements and global model predictions showed impacts 
on mesospheric HO x , NO x , and ozone. MLS observations of 
H0 2 displayed enhancements up to 1 ppbv due to the January 
2012 SPE period and up to 0.8 ppbv due to the March 2012 
SPE period in the northern polar mesosphere. The southern 
polar mesosphere was impacted as well during these SPE pe- 
riods, however, the increased solar radiation (and enhanced 
HO x loss) impeded the H0 2 increase in January. The SPE- 
caused H0 2 enhancements were fairly well simulated with 
the GSFC 2-D model. Both MLS measurements and GSFC 
2-D model predictions indicated middle mesospheric ozone 
decreases of > 20 % for several days in the northern polar 
region with maximum depletions > 60 % as a result of both 
the January and March 2012 SPE periods. Less SPE-caused 
ozone depletion was observed and simulated in the southern 
polar region in January 2012 due to the SPEs. Most of this 
mesospheric ozone depletion was the result of HO x catalytic 
loss cycles. 

Both ACE and M1PAS measurements showed large SPE- 
caused NO x enhancements (> 20 ppbv) in the northern polar 
mesosphere for about three weeks past the January 2012 SPE 
period. The southern polar mesosphere MIPAS observations 
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indicated a smaller SPE-caused NO x increase (<20ppbv) as 
a result of this January 2012 SPE period. The March 2012 
SPE period caused NO x measured enhancements > 20 ppbv 
and >10 ppbv for only several days in the northern and 
southern polar mesosphere, respectively. The GSFC 2-D 
model and the GMI 3-D CTM predicted NO x increases 
were slightly smaller (larger) than measured in the northern 
(southern) polar mesosphere as a result of the January 2012 
SPE period. The GSFC 2-D model and GMI 3-D CTM pre- 
dicted NO x enhancements as a result of the March 2012 SPE 
period were larger than observed in both polar hemispheres. 

The GSFC 2-D model and GMI 3-D CTM were used to 
predict medium-term impacts (~ months) from the January 
and March 2012 SPE periods. Stratospheric ozone was im- 
pacted in the polar regions by these solar events due to the 
transport of SPE-produced NO y to lower altitudes. The GMI 
3-D CTM predicted significantly larger upper stratospheric 
ozone destruction in the NH compared to the GSFC 2-D 
model, with a maximum of 13 %. In the SH, the model pre- 
dictions were generally similar, with a maximum ozone de- 
struction of 8-1 1 %. Total ozone reductions were predicted 
to be a maximum of 1—1.5 % in the SH polar region due to 
the SPEs and would not likely be observed in total ozone 
satellite or ground measurements. 
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